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Abstract. 2D hydrodynamical simulations are performed to examine the evap- 
oration and condensation processes of giant molecular clouds in the hot phase of 
the interstellar medium. The evolution of cold and dense clouds (T = 1000 K, 
uh = 3 cm~^, Af = 6 - 10^ Mq) is calculated in the subsonic stream of a hot tenuous 
plasma (T = 5 • 10^ K, uh = 6-10 cm ). Our code includes self-gravity, heating 
and cooling processes and heat conduction by electrons. The thermal conductivity 
of a fully ionized hydrogen plasma (k cx T^'''^) is applied as well as a saturated heat 
flux in regions where the mean free path of the electrons is large compared to the 
temperature scaleheight. 

Signiflcant differences occur between simulations with and without heat conduc- 
tion. In the simulations without heat conduction, the clouds outermost regions is 
stired up by Kelvin-Helmholtz (Klf) instability after only a few dynamical times. 
This prevents an infiltration of a significant amount of hot gas into the cloud before 
its destruction. In contrast, models including heat conduction evolve less violently. 
The boundary of the cloud remains nearly unsusceptible to KH instabilities. In this 
scenario it is possible to mix the formerly hot streaming gas very effectively with 
the cloud material. 
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1. Introduction 

The Interstellar Medium (ISM) can be described as a clumpy, inhomo- 
geneous three-phase medium (McKee & Ostriker, 1977: M077). The 
cold neutral phase (T ~ 80 K, n ~ 40 cm~^) is represented by the cores 
of molecular clouds which are enveloped by a warm neutral to slightly 
ionized medium (T ~ 8000 K, n ~ 0.3 cm~^). These two components 
are in pressure equilibrium if the gas is externally heated and can cool 
radiatively (Field et al., 1969). According to M077 they are embedded 
in a third hot, dilute phase of the ISM (T ~ 10*^ K, n ~ 10"^ cm'^). 
This component, which is produced by supernova (SN) explosions, can 
originally not be in pressure equilibrium with the colder phases and has 
therefore to expand vehemently. During this expansion, shocks arise, 
penetrate through and interact with the ambient clumpy ISM. In OB 
associations the sequential SNe accumulate and form a simple giant hot 
gas bubble, called superbubble, that expands preferentially towards the 
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halo and occupies a large volume (Norman & Ikeuchi, 1989). In order to 
allow for a more efficient cooling of the superbubble gas the evaporation 
of incorporated clouds are invoked (Silich ct al., 1994). This evaporation 
of cloud material should be caused by heat conduction of electrons from 
the hot to the cold gas. Our aim is to investigate the interaction between 
the streaming hot plasma of the ISM and the molecular clouds in order 
to determine the role of heat conduction during the cloud evolution. 

A situation described above comes into play in many astrophysical 
phenomena. As an example High- Velocity Clouds (HVCs) are radially 
falling towards the galactic disk (see Wakker & van Wocrden, 1997, for 
a recent review) and penetrate through the hot halo gas (Kerp et al., 
1996). Single HVCs are grouped together in complexes which can be 
found all over the sky. Interferometer measurements of the 21 cm Hi line 
at 1' resolution reveal the substructure of single complexes that consist 
of several small clumps embedded in larger emission regions (Wakker 
k. Schwarz, 1991). This, along with their line widths led Wakker & 
Schwarz conclude that the HVCs have a multi-phase structure consist- 
ing of a cold, dense core and a warmer, more tenuous envelope. Distance 
measurements for the cloud complexes remain difficult. For at least two 
of them upper limits for their distance are available (Danly et al. 1993; 
Keenan et al. 1995; van Woerden et al. 1997). These estimates reveal 
that they belong to the hot galactic halo. Although the origin of most of 
the complexes is still speculative (Blitz et al., 1999), HVCs are therefore 
a classic example for structures moving through a hot plasma. 

Another example for the scenario described above are the much more 
massive and larger speculative proto-globular cluster clouds (PCCs) 
with temperatures near 10'* K and densities several hundred times that 
of the surrounding gas. Therefore they are gravitationally unstable at 
masses larger than 10^ Mq . These clouds originate from condensations 
of thermally unstable gas with temperatures of some million Kelvin in 
the early epoch of galaxy formation and can be envisaged as progenitors 
of globular clusters (Fall & Rees, 1985). The PCCs had to resist their 
gravitational collapse for a sufficient time until star formation ignites. 
On the other hand, they are accelerated in the gravitational potential 
of the forming protogalaxy and therefore move through a hot plasma 
where they become subject to the growth of Kelvin-Helmholtz (KH) 
and Rayleigh- Taylor (RT) instabilities (Drazin &: Reid, 1981). 

Stable models for both, HVCs and PCCs, assuming hydrostatic and 
thermal equillibrium consist of large temperature and density gradients 
at the surfaces of the clouds where the energy densities of the hot, 
tenuous ISM and the warm, dense cloud become equal. Because of this 
large temperature gradient and the high temperature of the ambient 
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medium of the order of some million Kelvin, heat conduction has to 
play an substantial role on the evolution of such clouds. 



2. Heat conduction 



Heat conduction results in an additional heat flux g, that can be written 
in the classical case where the mean free path of the electrons is small 
compared to the temperature scaleheight as (Spitzer, 1962) 



with 
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where the Coulomb logarithm is 



InA = 29.7 + In 



lO^K 



(1) 
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with the electron density rig and the electron temperature T^, respec- 
tively. The diffusion approximation for the heat flux breaks down if the 
mean free path becomes comparable or even greater than the tempera- 
ture scaleheight. A common approach is to use a flux limited form, 
socalled "saturated" heat flux. This takes charge conservation into 
account and yields results in good agreement with more sophisticated 
treatments (e.g. Max et al., 1980) and with numerical simulations of 
laser heated plasmas (Morse & Nielsen 1973; Manheimer & Klein 1975). 
This saturated heat flux takes the form (Cowie & McKee, 1977: CM77) 



5 $s pc-^ 



(4) 



where c is the sound velocity and $s is a number less than or of the 
order of unity which embodies some uncertainties concerning the flux 
limited treatment (<&s = 1 in our calculation). In order to get a smooth 
transition between both the classical and the saturated regime we used 
the analytical form by Slavin & Cox (1992) 



(5) 



This guarantees that the smaller flux is taken if both differ significantly. 
As a criterion to separate both cases for a cloud of radius R embedded 
in a hot gas with temperature Tf, electron density ng/ and thus the 
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conductivity Kf, CM77 introduced a global saturation parameter which 
is essentially the ratio of the electron mean free path to the cloud radius 

0.08 KfTf 1.23 -lO^T^ 
(70 = = ^ . (6) 

$s Pfcj R Uef R 

For (To < 0.027/ $5 material condenses onto the cloud because radiative 
losses exceed the conductive heat input. For 0.027/$s < cro < 1 the 
clouds suffer classical evaporation, while for cjo > 1 the evaporation is 
saturated. The classical evaporation rate is given by CM77: 

while for the saturated case, rh-dass is multiplied by a function w{aQ) 
(Dalton &; Balbus, 1993) which lowers the mass loss rate for ctq > 0.01. 



3. Hydrodynamic Simulations 

The evolution of clouds in the subsonic stream of a hot plasma is studied 
by two-dimensional hydro dynamical simulations. For this the Eulerian 
equations are solved on a rectangular cylindrical symmetric "staggered 
grid" which is second order in space. The boundary conditions are 
chosen in a manner that material can flow out of the computational 
domain through the upper and right boundary. The physical parame- 
ters of the lower boundary that is also the symmetry axis are mirrored. 
The parameters, density and temperature, at the left boundary keep 
their initial values. Also the flow of the plasma is initialized by setting 
the velocity of the left boundary to the requested value. In order to 
follow the condensation of the streaming material onto the cloud a new 
quantity "colour" was introduced. The value of this new quantity is 
set in each cell to the density of hot ISM. At the beginning of the 
calculation only the cells not belonging to the cloud have a colour 
unequal to zero with a value equal to its density. During the calculation 
this quantity is advected like the others as e.g. density or energy density. 
For the advcction we used the monotonic transport of van Leer (1977). 
The Poisson equation for self-gravity was solved if the new evaluated 
gravitational potential differs from the old one by more than 10"^. 
The energy equation is extended by heating, cooling and heat conduc- 
tion, the latter one being the divergence of the heat flux (equ. 5). The 
adopted cooling function assumed collisional ionisation equilibrium and 
is a combination of the function introduced by Bohringcr & Hensler 
(1989) for T > lO^K and by Dalgarno & McCray (1972) for the lower 
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temperature regime. The heating function considers heating by cosmic 
rays, X-rays and the photoelectric effect on dust grains. 

3.1. Theoretical Model 

The initial profile of the cloud is generated for hydrostatic and thermal 
equilibrium. After setting the temperature Tism and particle density 
fiisM of the hot and tenuous outer medium the energy density Bism of 
the plasma is known. The density and temperature profile of the cloud 
is calculated by integrating the equations of hydrostatic and thermal 
equilibrium from inside out using the core temperature of the cloud 
as the inner boundary condition and is truncated when the energy 
density reaches ejsM- The model parameters are given in Table I. The 
parameters ctq and the evaporation timescale r^va are only defined for 
the simulation with heat conduction. The hot gas streams with 0.3 
Mach. The resulting (Jq for the simulation with heat conduction achieves 
a moderate saturation; the evaporation timescale Tcva is calculated from 
equation (7) correspondingly. The dynamical timescale as defined to be 
the sound travel time over one cloud radius, is 16.2 Myr. For the grid 
440 X 160 cells are used in the simulation with heat conduction and 
640 X 200 cells for the simulation without heat conduction. Thus, the 
spatial resolution amounts to 1.25 pc/cell. 



Table I. Model parameters of the simulation 



Tism 


nisM 


VlSM 


Rcld 




density 




Teva = M/rh 


(K) 


(cm--^) 


(km s^^) 


(pc) 


Mo 


contrast 




Myr 


5.6 • 10^ 


6.6 • 10"'' 


107.4 


41.3 


6.4 • 10* 


1.2 • 10* 


10.5 


279 



4. Results 

The inclusion of heat conduction has the tendency to stabilize massive 
and large clouds. This is illustrated comparing the density contours at 
the same times (30 Myr and 50 Myr after the beginning of the calcu- 
lations) for the simulations without (Fig. 1) and with heat conduction 
(Fig. 2). Without heat conduction the edge of the cloud is torn by KH 
instability already shortly after the beginning when the density contrast 
is lowered. Also a very complex velocity structure with many vortices 
behind the cloud has formed. On the other hand, heat conduction sup- 
presses large scale KH instabilities and only a single large circulation in 
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the slipstream of the cloud is visible. This behaviour can be understood 
if the conducting electrons move along the temperature gradient, i.e. 
radially toward the cloud. This leads to a deceleration of the contact 
interphase so that the relative velocity decreases below a critical value 
that determines the onset of KH instabilities. Also the inner parts of the 
clouds develop differently. While in both simulations the cores remain 
very dense, the model without heat conduction develops a compound 
envelope, with a radially decreasing density and a diluted outer part of 
material that is only slightly gravitationally bound. 



181.8 km/s 




200 300 pc 400 500 

Figure 1. Density contours for the simulation without heat conduction at 30 Myr 
(upper panel) and 50 Myr (lower) after the beginning of the calculation. The con- 
tour lines represent 5, 10, 50, 100,. . . xpisM- The velocity arrows scale linearly with 
respect to the maximum velocity shown in the upper right. 
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Figure 2. Density distribution for the simulation with heat conduction at 30 Myr 
(upper panel) and 50 Myr (lower) after the beginning of the calculation.The contour 
levels are the same as in Fig. 1 



The cloud model with heat conduction remains more compact with 
bound edges. In regions where the heat conduction exceeds the cooling 
the density distribution is homogenized while in very dense regions like 
e.g. near the core the density structure is unaffected by heat conduction. 
The tendency to smooth out temperature and density inhomogenities 
is also the reason why KH instability is suppressed. At the edge of 
the cloud an interface builds up which broadens the transition region 
from the cloud to the outer medium. This transition region is visible 
in Fig. 2 as a broad band of contour lines where the former steep 
temperature and density gradients are diminished. Within this interface 
at the cloud's surface the sheer velocity is therefore lowered and the KH 
instability is reduced, by this, enabling material from the flow to settle 
onto the cloud. This difference is clearly discernible in Fig. 3. After 76 
Myr almost 160 Mq have stuck into the cloud what amounts to almost 
0.3% of the total cloud mass. The simulation without heat conduction 
shows an accretion that is only 1/3 of that one with heat conduction 
because the accreted material is stripped off from the cloud again before 
it can be mixed with deeper layers of the cloud. The evolution of the 
total mass (see Fig. 4) differs only marginally. Long-term simulations 
have to be performed in order to investigate whether the total mass 
reaches a constant value. Nevertheless, it is obvious, that the cloud 
evolution is not suffering evaporation at that stage as claimed by CM77. 
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The mass-loss rate corresponding to the analytical expression given by 
equ. 7 is by far too high. The question has to be addressed whether the 
assumptions by CM77 arc vahd in such a more realistic case in which 
instabilities occur and a streaming hot plasma serves as an infinite 
reservoir. In order to answer it, a parameter study covering both the 
condensation and evaporation regime proposed by CM77 has to be 
performed and investigations must include a higher spatial resolution. 



cn 
o 


0.0030 




without heat conduction/, 

with heat conduction «/ 




E 










o 










o 


0.0020 








J mass 


0.0010 








0) 




















acc\ 


0.0000 









20 40 60 

Mvr 



Figure 3. Evolution of the accretion of the streaming material onto the cloud. 
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Figure 4- Evolution of the gravitationally bounded mass with the regard to the 
initial mass. 
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